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EVALUATION RESULTS AND ACCURACY IMPROVEMENT 
OF TWO FLUX-SWITCHING MATHEMATICAL MODELS 
FOR MAGNETIC CORES 

By J. Larry Spencer 
Langley Research Center 

SUMMARY 

The results of an evaluation of two flux-switching models for magnetic cores, the 
Nitzan model and the Betts and Bishop model, are presented. The two models are com- 
pared on the basis of the types of core switching each describes and the fidelity with which 
each predicts the output voltage waveform of a ferrite core with inside diameter of 50 mils 
(1.27 mm) and outside diameter of 80 mils (2.03 mm). The results indicate that the 
Nitzan model is more accurate and more nearly complete because it describes shuttling 
and partial switching. The Nitzan model is more applicable to cores of the same material 
but different sizes because it uses separate parameters related to the core materials and 
the core dimensions. A method is presented by which both models can be extended to 
describe more accurately the switching of cores with smaller ratios of inside to outside 
diameter. 


INTRODUCTION 

For the mathematical analysis and design of magnetic core circuits, it is necessary 
to have a mathematical model to describe the operation of a magnetic core. A flux- 
switching model for a magnetic core describes mathematically the output voltage-time 
waveforms appearing on the sense winding as a function of the current through the drive 
winding. A core model for use in circuit analysis and design should accurately predict 
the output voltage of a core under all switching conditions but be simple enough to aid the 
designer of magnetic circuits. 

Most of the existing core models are based on an approach developed by Menyuk 
and Goodenough (ref. 1), which predicts the output voltage on the sense winding as a func- 
tion of the input current and initial flux state of the core. Seven flux- switching models 
for magnetic cores were evaluated by Hesterman (ref. 2), who concluded that the linear 
parabolic flux-switching model is the best for fast switching and the nonlinear parabolic 
model is the best for slow switching. Nitzan et al. (refs. 3 to 10) combined the linear 
and nonlinear parabolic models and extended the model to describe shuttling and incomplete 



switching. In Hesterman's report, reference 2, a model developed by Betts and Bishop 
(ref. 11) was listed as promising, but the model was not evaluated. In the present study 
the models developed by Betts and Bishop and by Nitzan are compared and evaluated. 

These models describe switching of a magnetic core from the maximum or minimum 
residual-flux level. A method is also presented for improving the accuracy of both models 
in predicting the output voltage of a "thick" core. The symbols used are defined in appen- 
dix A and a brief description of the models is presented in appendix B. 

COMPARISON OF THE MODELS 

The Betts and Bishop model and the Nitzan model are of a complexity that dictates 
their solution on a computer. The computer times for the two models are nearly equal, 
and therefore the speed of solution is not a factor in their comparison. The two models 
are compared first on the basis of the types of core switching each describes and then on 
the basis of the accuracy with which each model predicts the output voltage of a core. A 
typical hysteresis loop for a magnetic core is shown in figure 1. The Betts and Bishop 
model describes only the complete switching (±4> r to of a magnetic core. The Nitzan 

model describes complete switching and also partial switching (±<f> r to £+$ s ) because 
the static hysteresis loop and the nonlinear region of the peak output voltage curve are 
included in its parameters. The output voltage produced when the core is shuttled (+<h r 
to ±$ s ) is Predicted by the elastic switching part of the Nitzan model. The Nitzan model 
is more nearly complete than the Betts and Bishop model. Neither model will describe 
switching from a partially set state - that is, switching from some flux state other than 
+<i> r . 

The accuracy of the models was studied by comparing the output voltage predicted 
by both models with the experimentally determined output voltage under the same drive 
conditions. Four different drive conditions were used. The ferrite core, which had an 
inside diameter of 50 mils (1.27 mm), an outside diameter of 80 mils (2.03 mm), and a 
height of 30 mils (0.76 mm), was completely switched by each drive condition. (This 
core is called the 50-80 core herein.) The test circuit shown in figure 2 was used to 
determine the experimental output voltage. The generator current waveform is shown in 
figure 3. The measurement of the experimental parameters for the models is described 
in detail in reference 12. Also in reference 12 is a complete set of computer programs 
in FORTRAN IV that describe the Nitzan model and the Betts and Bishop model. A 
CDC 6600 computer was used to obtain the predicted output voltage of the models for the 
accuracy comparison. 

The voltage predicted by the models was compared with the experimentally obtained 
voltage for the following driving currents (pulse 1 in fig. 3): 
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Case 

Current, 

^ S ’ 

amperes 

Rise time, 
trise* 
/isec 

Ramp current switching 

0.79 

2.95 

Ramp-step current switching 

.77 

.88 

Slow step current switching 

.79 

.04 

Fast step current switching 

1.6 

.035 


Pulse 2 is used to reset the core. In ramp current switching the core is switched by a 
ramp driving current, whereas in ramp-step current switching the core is partially 
switched by a ramp and the switching is completed by a constant current. In slow step 
and fast step current switching the core is switched by a constant current having a very 
short rise time. The output voltage waveforms predicted by both models and the experi- 
mental output voltage waveform are shown in figures 4, 5, 6, and 7 for ramp, ramp-step, 
slow step, and fast step switching, respectively. 

The Nitzan model describes the output voltage waveform for a core when partially 
switched. Figure 8 shows the waveform of the output voltage predicted by the Nitzan 
model and the experimental output of the 50-80 core when 11.5 maxwells (115 nano- 
webers) of the 12.6 maxwells (126 nanowebers) of flux is switched by a driving current of 
0.4 ampere with a rise time of 0.04 microsecond. The Betts and Bishop model will not 
predict partial switching. 

The Nitzan model also describes shuttling of the core. To apply the model to a 
shuttled core, the time rate of change of the applied current must be determined. The 
driving current was fitted by a polynomial as shown in figure 9. The equation of the 
driving current was differentiated with respect to time and multiplied by a core constant e 
to obtain the prediction of the shuttling output voltage by the Nitzan model. The accuracy 
of the prediction is shown in figure 10. 

EXTENSION OF THE MODELS 

In figures 4 to 7 the calculated output voltage curves from both models are shifted 
to the right and up from the curve of the measured voltage. The models are more inaccu- 
rate when describing switching by a ramp driving current than by a step driving current. 
The main source of inaccuracy of the models for these particular tests is the assumption 
that the 50-80 core has a diametral ratio close to unity (thin). Restricting the models to 
thin cores justifies the assumption that the magnetic field intensity throughout the core 
material is a constant. The actual magnetic field intensity in a magnetic core is inversely 
proportional to the distance from the center of the core. In reference 13 the Betts and 
Bishop model was extended for cores of any diametral ratio by determining the output 
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voltage of a thin core in terms of the radius and analytically integrating this voltage from 
the inner to the outer radius of the core. The integration of the Betts and Bishop model 
complicates the model so that it cannot be used practically in circuit analysis. 


By analytically dividing a thick core into concentric rings, where each ring can be 
assumed to be a thin core, and applying the thin-core model to each ring, the models are 
extended to describe switching of a thick core. The output voltage-time response of a 
thick core, as predicted by the extended models, is the analytical sum of the output 
voltage-time responses from each of the concentric rings. The two extended models, 
using two rings, were applied to the 50-80 core subject to the same test conditions as in 
figures 4 to 7. The results are shown in figures 11 to 14. (The scale changes on the 
time axis of the figures for the two-ring models result from the automatic plotting of the 
figures.) In all four cases, the accuracy of both models was improved. The greatest 
improvement exists for switching when the driving current is a ramp (compare figs. 4 
and 11). By comparing the accuracies with which the models predict ramp switching and 
fast step switching in figures 4 and 7, respectively, it can be seen that what is considered 
to be a thin core for fast switching is a thick core for slow switching. In figure 11 the 
first peak in the output voltage waveform is caused by the domination of the faster 
switching, higher amplitude output voltage of the inner ring and the second peak is caused 
by the slower switching, lower amplitude output voltage of the outer ring. If more than 
two rings are used, the difference in the output voltages from the rings is less apparent 
in the sum. Therefore, the total output voltage obtained by summing the predicted output 
voltages of the rings is a smoother, more accurate prediction of the experimental output 
voltage of the core. 

To prove that additional rings will improve the accuracy, a 10-ring Nitzan model 
was applied to ramp-current switching and partial switching, the two cases in which the 
one-ring model is most in error. Figure 15 shows the prediction of the output voltage 
for ramp-current switching with a 10-ring model. There is more improvement in the 
accuracy of the model in going from one to two rings than in going from two rings to 
10 rings. The rings also improve the accuracy of the Nitzan model when describing 
incomplete switching. The Nitzan model for describing partial switching with 10 rings 
is shown in figure 16. The output voltage of a shuttled core as predicted by the Nitzan 
model is not affected by dividing the core into rings. 


The percent of average deviation of the output voltage by the model from the experi- 
mental output voltage was calculated for both models under each driving current from the 
following equation: 


Percent average 


r W [W> - V °,m(t>]* 

deviation = — 


sw^o,p,exp 


x 100 
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Table I gives the percentage of average deviation from the experimental output voltage 
for each case of complete switching and partial switching. The data show that for a par- 
ticular core and a given drive condition, the accuracy of the models can be improved by 
using a larger number of rings to represent the thick core. However, the time required 
for computation increases in direct proportion to the number of rings used. Greater 
improvement is achieved with the first segmentation of the core into rings, and successive 
subdivisions result in less improvement. 

TABLE I.- COMPARISON OF RESULTS FROM MODELS AND EXPERIMENT 


Percentage of average deviation from the 
experimental output voltage for — 


Case 

Betts and Bishop 
model with — 

Nitzan model with — 


1 ring 

2 rings 

1 ring 

2 rings 

10 rings 

Ramp current switching 

44.8 

31.5 

22.5 

11.7 

n.i 

Ramp -step current switching 

22.1 

17.9 

12.1 

7.3 

Not tried 

Slow step current switching 

10.6 

8.6 

8.5 

5.4 

Not tried 

Fast step current switching 

6.9 

6.4 

7.1 

6.7 

Not tried 

Partial switching 

Not applicable 

88 

Not tried 

13.5 


CONCLUDING REMARKS 

An evaluation of two flux -switching models for magnetic cores, the Betts and Bishop 
model and the Nitzan model, has been made. It is concluded that the Nitzan model is a 
more nearly complete model because it will describe shuttling and partial switching, 
whereas the Betts and Bishop model will not. Neither model will describe switching from 
a partially set state, which is an important switching condition in magnetic circuits. An 
advantage of the Nitzan model is that it uses one set of parameters related to the material 
of which the core is constructed and another set related to the physical dimensions of the 
core, whereas the Betts and Bishop model uses parameters of these quantities combined. 

A method has been introduced for improving the accuracy of both models in pre- 
dicting the output voltage of a "thick'* core. It was found that the accuracy of both models 
in predicting the output voltage of a thick core is improved by accounting for the magnetic 
field intensity variation through the core. The extended models accounted for the field 
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variation by analytically segmenting the core into concentric rings and applying the model 
to each ring. The number of segmentations needed for accurate description of the output 
voltage is dependent upon the diametral ratio of the core, the driving current, and the 
squareness of the hysteresis loop of the core material. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., September 4, 1970. 
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APPENDIX A 


A 

B 

Bp 

Bp(H) 

B r 

B s 

Cn 

f(r) 

H 

Hl-2 

H 2 -3 

H a 

Hdyn 

H ext 

H n 


SYMBOLS 

cross-sectional area of the core 
magnetic flux density 

first derivative of peak flux density with respect to time 

equation relating the peak flux density to the driving magnetic field intensity 
maximum residual flux density 
saturation flux density 

coefficients of the characteristic equation for the Betts and Bishop model 
characteristic equation of normalized output voltage waveform 
driving magnetic field intensity 

magnetic field intensity boundary between region 1 and region 2 of the 
Bp(H) curve 

magnetic field intensity boundary between region 2 and region 3 of the 
Bp(H) curve 

magnetic field intensity asymptote for a parabola in describing the 
static B-H loop 

dynamic magnetic field intensity threshold 

magnetic field intensity obtained by linearly extrapolating the Bp(H) curve 
to the H-axis 

magnetic field intensity asymptote for a parabola in describing the static 
B-H loop 
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APPENDIX A — Continued 


Hq magnetic field intensity at which the parabola in describing the static B-H 

loop intersects the B = -B s line 

H s t static magnetic field intensity threshold 

h axial dimension of the core 

I amplitude of step driving current 

If _2 current boundary between region 1 and region 2 of the 4>p(I) curve 

^2-3 current boundary between region 2 and region 3 of the <lp(I) curve 

Idyn dynamic current threshold 

I ex t current threshold obtained by linearly extrapolating the third region of the 

peak output voltage curve 

Ii input driving current 

I ss steady-state driving current 

I s t static current threshold 

l magnetic path length 

f mean magnetic path length, 7r(ri + r Q ) 

l ^ inside circumference of the core 

l Q outside circumference of the core 

Rp slope of peak output voltage curve 

inside radius of the core 

r Q outside radius of the core 

t time 
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APPENDIX A — Continued 


tp time at which peak output voltage occurs 

t r ise rise time of the current 

t sw switching time 

tt time at which the driving current reaches the threshold of the core 

Vo, exp 00 experimental output voltage 

V Q) m output voltage predicted by a model 

V Q p peak output voltage 

Vo,p,exp peak experimental output voltage 

W c core constant equal to V 0> p >ex ptp 

6 coefficient of inelastic B 

e coefficient of elastic switching 

k i coefficient of inelastic B in nonlinear Bp(H) for Hgj. = H < Hj_2 

«2 coefficient of inelastic B in nonlinear Bp(H) for Hj_2 = H < H2_3 

Xj coefficient proportional to <i>p(I) for I st ^ I < Ij g 

^2 coefficient proportional to 4 >p(I) for Ii_2=I < l2-3 

exponent of 4 > p (I) for I st S I < I^_ 2 
i^2 exponent of ip(I) for Ij_2 = I < I2-3 

p coefficient of inelastic 4 

t normalized time, t/tp 

$ magnetic flux 
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APPENDIX A - Concluded 


$ 

*el 

^inel 

^max 

*p(I) 

*s 

♦st® 


first derivative of flux with respect to time (voltage) 
elastic portion of output voltage per turn of a core 
inelastic portion of output voltage per turn of a core 
maximum flux switched, 3> s + <f> r 

equation relating the peak output voltage to the step driving current 
maximum residual flux 
saturation flux 

equation of the static hysteresis loop 
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APPENDIX B 


DESCRIPTION OF MODELS 


The Betts and Bishop Model 

The Betts and Bishop model mathematically predicts the output voltage of a thin 
square-loop magnetic core which is excited by a drive current. The Betts and Bishop 
model describes only complete switching of a core; that is, switching from +<l> r to +$ s , 
which means 

Vo.expW dt = *max = Constant (1) 

The Betts and Bishop model is based on two assumed conditions: (1) The open-circuit 

voltage-time waveforms caused by step input currents are identical when normalized with 
respect to peak amplitude and the time at which this occurs, and (2) the peak output volt- 
age of an unloaded core is directly proportional to the magnitude of the step input current 
above the current threshold of the core I e xt* The normalized output voltage waveform 
is known as the characteristic curve, and the equation of this curve is the characteristic 
equation f(r). The equation of the peak output voltage curve is 

v 0 ,p = Rp(I - ^ext) (2) 

and the characteristic equation is 
N 

f(T) = £ C n T n (3) 

n=0 


The descriptive equations of the Betts and Bishop model for a step input current are 


Vo,exp(t) ~ 


T = 


[%>(I “ *ext)f( T ) 



" Iext) I 


(l < Iext) 

(I ^ Iext) 

(l < Iext) ^ 
(l= Iext) 


(4) 


where W c is the product of a peak output voltage and the time at which the peak 
occurred. 

The output voltage predicted by the model is the amplitude of the characteristic 
equation scaled according to the peak output voltage curve, with time scaled according to 
the driving current and flux capacity of the core. The descriptive equations for the Betts 
and Bishop model are extended to describe switching with an arbitrary driving current by 
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APPENDIX B - Continued 


using a series of step driving currents to represent the arbitrary driving current and 
letting the time duration of each step approach zero. The Betts and Bishop model for an 
arbitrary driving current is 


^o,exp^) 



(Ii(t)<Iextf 
(Ii(t) = I e xt) 



(li(t) < Igxt) 

Ijp- f - I ext (t - t t ) (Ii(t) = I ext ) 

c tj- 




(5) 


The time at which the driving current reaches the threshold of the core is denoted by t^. 


The Nitzan Model 

The Nitzan model predicts the output voltage of a thin core when completely switched, 
incompletely switched, and shuttled. A detailed description of the Nitzan model, which is 
only summarized in this report, is in reference 5. The model predicts the output voltage 
as the sum of the elastic switching component and the inelastic switching component: 


^o,m - $ e l 


+ 


inel 


The elastic flux switched is the flux that will return to its initial value when the applied 
current is reduced to zero. The inelastic flux switched is the part of the total flux that 
will not return to its initial value when the driving current is reduced to zero. Inelastic 
switching produces the square hysteresis loop of a magnetic core. 


The Nitzan model for elastic switching is 
dl (t) 


The slope of the static hysteresis loop in the saturated region is denoted by e 
approximated by 


( 6 ) 

and is 


e 





(7) 


Only the elastic switching part of the Nitzan model is used to predict shuttling of the core. 

The Nitzan model for inelastic switching is defined by the nonlinear differential 
equation 
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APPENDIX B - Continued 


•taelW = *pW “ 


2$(t) + <& s - •rtOT 


$s + *st(0 


(8) 


where 4>p(I) is the equation for peak output voltage as a function of amplitude of the step 
driving current, and «& s t(I) is the equation describing the static hysteresis loop of the 
core. The static hysteresis loop is the relationship between current and flux of a mag- 
netic core when switching as slowly as possible. Equation (8) was semiempirically 
derived by fitting a parabola to the relationship between voltage and flux of a magnetic 
core. 


The equation of the static hysteresis loop $ s t(I) is derived in such a way that it 
is represented by a set of parameters dependent on the core material and a set dependent 
only on the physical dimensions of the core. The advantage of having these two inde- 
pendent sets of parameters is that after experimentally testing one core and calculating 
the material parameters, the model can be applied to any core of this same material by 
merely adjusting the physical-dimension parameters in the model. The loop of flux 
density as a function of magnetic field intensity (B-H loop), from which the <h-I loop 
is derived, depends only on the magnetic material. The B-H loop of the core being 
modeled is analytically described by hyperbolas as in figure Bl. By means of the 
relationships 

$ = J B • dA 

and 

I = (j) H ’ d£ 

the equation describing the B-H loop for a material is converted to equations for the 
$-1 loop for a core constructed of that material. The equations of the static $-1 loop 
for a toroidal core are 


*s-*r .... 

(l 0 - ?l)H a I - H a /j 


( ! s I*) 


(*s + *r) H q 

I 

(^o “ ^i) H n 

Hq 

(*s + ^r)Hq 

7 

(^o “ ^i)Hn 

L o 


1 _ _ 1 _ 
H n "H q 


1 - 5 * 

H q 

H n ii 


7 7 . t ( 1 1 \ 1 ~ H n^o 

Z °-^ +I (H n 'H q j ln i-H n Z i 


- $ r (I st < I < HqZ 0 ) ) (9) 




(Hqi 0 = i) 
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APPENDIX B — Continued 



Figure B2.- Peak output voltage curve with regions. 


The relationship of peak output voltage to amplitude of step driving current is shown 
in figure B2. The curve is divided into two nonlinear regions and one linear region. The 
equations describing the regions are as follows: 


*p(I) = < 


f' 

0 

(0 i 1 < 1st)' 

Hi 1 ~ 1st) 11 

(1st = I < I 1 - 2 ) 

*2(1 - Idynf 2 

(I 1 - 2 SKI 2 - 3 ) [ 

p(I " lext) 
v_ 

( : 2-3 s >), 


( 10 ) 


A computer program fits equations (10) to the experimental peak output voltage data 
by adjusting Xj, X 2 , v^, i/g, p, I ex t, Idyn> 311(1 1st- If the core is assumed to be a 
thin toroid, the material parameters of the core are determined by use of the following 
equations: 
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APPENDIX B - Concluded 


H, 


dyn 




H 


x 2-3 
2-3 = — 


6 = 


B s = 


l 

h(r Q - H) P 

*s 

h ( r o - r i) 


H ex t = ^ 

l 


h *‘-T 


H 


1-2 


l l-2 


B r = 




h ( r o " r i) 


f 2 


K 1 = TJZ TT a 1 

1 h ( r o - r i) 


The parameters H dyn , H 2 _ 3 , H ext , H st , H 1 _ 2 , 6, B s , B r , k 2 , and Kj are 

material parameters. 

Because of the complexity of the peak output voltage equations (10) and static hys- 
teresis loop equations (9), the descriptive equation (8) must be solved numerically on a 
digital computer. An arbitrary driving current is represented by a series of step driving 
currents of very small time duration, and a simple predictor-corrector method is used 
to solve equation (8) for •f'inelW- By summing $ inel and $ e j, the total output voltage 
of a core driven by an arbitrary drive current is predicted by the Nitzan model. 
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Figure 2.- Test circuit. 



Figure 3 .- Pulse sequence from the pulse generator of figure 2. 
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Figure 6.- Output voltage for slow step switching. 
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e ‘T he aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and Notes, 
and Technology Surveys. 
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NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



